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Nonclassical states of light are important resources 
for quantum metrology ^ [2] , secure communication 3J , 
quantum networks [IHSj, and quantum information pro- 
cessing [71 |B] . Of particular interest for these applications 
are traveling wavepackets prepared with a definite num- 
ber of photons in a continuous temporal mode, known as 
continuous-mode Fock states [9Ul2j. As the generation 
of such states becomes technologically feasible p!5H^ a 
theoretical description of the light-matter interaction [2S] 
becomes essential, see Fig. [T] 

Previously, aspects of continuous-mode single-photon 
states interacting with a two-level atom have been ex- 
amined. Others have investigated master equations 
two-time correlation functions [SSI El]; properties of scat- 
tered light [27t;35. ; and optimal pulse shaping for excita- 
tion [55H55] . The results in these studies were produced 
with a variety of methods which have not been applied to 
many systems other than two-level atoms or Fock states 
where ^ 1. 

One way to approach such problems is through the 
input-output formalism of Gardiner and Collett [40Vl44j . 
A central result of input-output theory is the Heisenberg- 
Langevin equation of motion driven by quantum noise 
that originates from the continuum of harmonic oscilla- 
tor field modes [45] . The application of input-output 
theory to open quantum systems has historically been 
restricted to Gaussian fields [U HBj — vacuum, co- 
herent, thermal, and squeezed — with several notable ex- 
ceptions [26, 47 -50^. 

In this article we present a unifying method, based on 
input-output theory, for describing the interaction be- 
tween a quantum system and a continuous-mode Fock 
state. Consequently our formalism encapsulates and 
extends previous results. Specifically our method al- 
lows one to derive the master equations and output 
field quantities for an arbitrary quantum system interact- 
ing with any combination of continuous-mode A^-photon 
Fock states. 

This article is organized as follows: In Sec. [l] we intro- 
duce the white noise Langevin equations of motion, the 
mathematical description of quantum white noise, and 
the formal definition of continuous-mode Fock states. In 
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FIG. 1: Schematic depiction of a traveling wavepacket in- 
teracting with an arbitrary quantum system. The tempo- 
ral wavepacket is described by a slowly- varying envelope ^(t) 
which modulates fast oscillations at the carrier frequency. We 
consider the case where the wavepacket is prepared in a non- 
classical state of definite photon number. 



Sec. [IT] we present the first main results: the method for 
deriving master equations for systems interacting with 
continuous-mode Fock states and related output field 
equations. This result is then extended in Sec. |III| to 
continuous-mode "A-photon states," where the spectral 
density function is not factorizable. Then, in Sec. IV we 



apply our formalism to the study of a two-level atom in- 
teracting with wavepackets prepared in A^-photon Fock 
states. This application is intended to serve as an in- 
structive example that reproduces and extends results in 
previous studies [36H38] . In Sec.|v] we present the second 
main result: master equations and output field quanti- 
ties for a system interacting with Fock state wavepackets 
in two modes (e.g. spatial or polarization). This sets 
the stage for the study of many canonical problems in 
quantum optics. As a two-mode example, we examine 
the scattering of Fock states from a two-level atom in 



Sec. VI Finally, we conclude in Sec. VII with discussion 
and possible applications. 



I. MODEL AND METHODS 

A description of a system interacting with a travel- 
ing wavepacket naturally calls for a formulation in the 
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time domain. The input-output theory developed in the 
quantum optics community provides such a description 
[HHm liSHlS] . Often input-ouput theory is formulated 
for a one-dimensional electromagnetic field, although this 
is not a necessary restriction [35]. (Such effective one- 
dimensional models are typically thought about in the 
context of optical cavities [511 or photonic waveguides 
[35, 52 54J.) In this formalism the rotating wave approx- 
imation, the weak-coupling limit (the Born approxima- 
tion) , and the Markov approximation are made |551 156j . 
Strict enforcement of these approximations is known as 
the quantum white noise limit |57j . 

In Appendix ( A 1 ) we review the quantum white 



noise limit; other introductory material can be found in 
Refs. im [421 ISSl [58]. The main result is a quantum 
stochastic differential equation (QSDE) for the unitary 
time evolution operator that governs the system-field dy- 
namics. From this equation one can derive QSDEs for 
system and field operators driven by white noise, also 
known as white noise Langevin equations. It is these 
equations of motion that lie at the heart of the deriva- 
tion of Fock-state master equations. 

The Langevin equations derived in the white noise 
limit are in Stratonovich form [9l [42l [59] . Stratonovich 
QSDEs obey the standard rules of calculus, but expecta- 
tions can be hard to calculate because the quantum noises 
do not commute with the operators to which they couple. 
Stratonovich QSDEs can be converted to an equivalent 
form known as the Ito QSDEs. In Ito form the quan- 
tum noises commute with the operators to which they 
couple, which facilitates taking expectations. However, 
differentials must be calculated to second order |42]. To 
derive master equations we will be taking expectations 
over field states and consequently will work solely with 
Ito QSDEs. 



A. Derivation of the vacuum master equation from 
the Ito Langevin equations 

Consider an arbitrary system operator in the interac- 
tion picture, X{t), with the initial condition X{to) — 
X /field- The time evolution of X is given by the Ito 
Langevin equation [see Appendix (A3|] 



dX ^{i[H,X]+ C^L]X)dt 

+ [L^ , X] SdBt + [A, L\dBl 
+ {S^XS ~ X)dAt, 



(1) 



where the action of the superoperator is 

C^L]X = L^XL- ^{L^LX + XL^L) . (2) 

The operators (5, L, H) act on the system Hilbert space. 
The quantum noise increments dBt, dBt\ and dAf are 
field operators, discussed in more detail shortly. 

The first two terms in Eq. ([T]) describe smooth evo- 
lution from an external Hamiltonian on the system and 



from a Lindblad-type dissipator. The second two terms 
describe the infiuence of quantum noise through coupling 
of a system operator L linearly to the field operators, 
e.g. dipole-type coupling. The final term arises from 
coupling of a system operator to a quantity quadratic 
in the field operators, such as photon number. Such ef- 
fective couplings appear in optomechanical systems |60| 
and arise after adiabatic elimination of the excited states 
in multi-level atoms 61 , for example. 

Let us return to the discussion about the quantum 
noise increments di?t, dBt\ and dAf. These field opera- 
tors are defined in terms of the fundamental field opera- 
tors h{t) and W (t) whose time arguments are mode labels 
rather than indicators of time evolution. They are often 
referred to as white noise operators because they satisfy 
the singular commutation relations [6(s), b^{t)] = S{t — s). 
This is akin to classical white noise which is (5-correlated 
in time. Due to the singular nature of b(t) and b\t), it is 
preferable to work with the quantum noise increments: 



/t+dt /'Z+az 
dsb{s) and dBl= dsb\s), (3) 

(4) 



pt+dt 



It Jt 

pt+dt 

dAt = / dsb\s)b{s) 



which drive the Heisenberg dynamics in Eq. ([T]). 

Under vacuum expectation, the calculus rules for ma- 
nipulating QSDEs are summarized by the relations 



dBtdBl = dt, dBtdkt = dBt, 
dktdkt = rfAt, dktdBl = dB^. 



(5) 



These composition rules are often referred to as the vac- 
uum Ito table. 

As a prelude to the derivation of the Fock-state master 
equations, we derive the vacuum master equation. First, 
we take vacuum expectations of Eq. p]) u sing the fol- 
lowing notation (to be explained in Sec. [II[): Eq^oMA] — 
Tr[(psys«) |0)(0|)^dA]. Consequently, we need the action 
of the quantum noise increments on vacuum, 



dBt\Q) =0, 
dAtiO) =0. 



(6) 
(7) 



All of the quantum noise terms in Eq. ([T]) vanish under 
vacuum. Then, using the cyclic property of the trace we 
obtain the vacuum master equation: 



dt 



(8) 



where the Lindblad superoperator is defined as 

£[L]q = LgL^ - l{L^Lg+ gL^L), (9) 

and the subscripts on gQ Q denote that Eq. ([8| is a vacuum 
master equation. 
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B. Continuous- mode Fock states 

A continuous-mode single-photon state O [TOl [12] can 
be interpreted as a single photon coherently superposed 
over many spectral modes |62l 163] with weighting given 
by the spectral density function (SDF) ^(w), 



/ dc.|H6tH|o). 



(10) 



We focus on quasi-monochromatic wavepackets, where 
the spectral spread is much smaller than the carrier fre- 
quency, Ao; ^ u)c [83J. This holds for optical carriers, 
whose bandwidths are small relative to the carrier fre- 
quency. Then, we can define a slowly-varying envelope 
^(cj) rotating at the carrier frequency. 



(11) 



where lOc is near any relevant system frequencies. 
The Fourier transform of the slowly-varying envelope, 
= ^{t), characterizes a square-normalized tem- 
poral wavepacket, J dt\S^{t)\'^ = 1. In the time domain, 
and within the quasi-monochromatic approximation, the 
single-photon state in Eq. (ITo]) becomes [TUj . 



\k) = J dsaspim 



(12) 

where we have absorbed the possible detuning from the 
system frequency into £,{t). The operator B^{£_) creates a 
single photon in the wavepacket ^(i). Equation ( [T2| ) can 
be interpreted as a superposition of instantaneous pho- 
ton creation times weighted by the temporal wavepacket. 
Since the white noise operators are defined in the inter- 
action picture, it is clear that ^(t) is a slowly-varying 
temporal envelope rotating at the carrier frequency. By 
focusing on quasi-monochromatic wavepackets we ensure 
the approximations made in the quantum white noise 
limit are not violated. 

A straightforward extension leads to the definition 
of normalized, continuous-mode Fock states (referred to 
hereafter as Fock states) in the wavepacket ^(i) with N 
photons pjj. 



1 



ds^ispis 



N 



|0) 



(13a) 
(13b) 



The Fock states in Eq. ( |13[ ) are a subset of more general 
A^-photon states for which the SDF is not factorizable 
[n]. In Sec. 



state. The derivation is performed in the interaction 
picture where the time-dependent operators evolve ac- 
cording to Eq. ([1]). To facilitate the derivation we first 
introduce notation convenient for representing expecta- 
tions with respect to a particular field state. It should be 
noted that our method is a generalization to A^-photon 
states of a method introduced in Refs. [321110] for a single 
photon. 

Assuming no correlations before the interaction, the 
total system is described by the product state 



pito)=Psys(^\N^){N^l 



(14) 



with the system in the state psys and the field in the Fock 
state \N^)- Using the Hilbert-Schmidt inner product for 
operators A and B, 



{A\B) = Ti[A^B], 



(15) 



one can take expectations with respect to system and/or 
field states. For the following derivation it is necessary 
to define the asymmetric expectation value, 



IEm,n[0] = Trgys+ficid (Psys |n^)(m^|)^ O 



(16) 



where O is a joint operator on the system and field 
and is not necessarily separable. We use a convention 
where capital letters, |A^^) denote the number of pho- 
tons in the input field. Lowercase letters, e.g. where 
n = {0, A^}, label "reference" Fock states to which the 
system couples. Using the Hilbert-Schmidt inner prod- 
uct, we define a set of generalized density operators Pm,„, 
first introduced in Ref . |26j , by tracing over only the field 



in Eq. (16) 



(17) 



Such generalized density operators were also used in 
Refs. [311[S0] for a single photon. We delay the interpreta- 
tion of these generalized density operators until Sec. |II A| 
As the trace in Eq. ( 16 ) is over both system and field. 



it gives a c-number expectation value. Using the partial 
trace we also define an asymmetric partial expectation 
over the field alone which results in an operator. We 
define this operation with the notation [Mj, 



(O) = Trfioid 



(/sys® |n5)(TO^|)^0 



(18) 



We base our derivation on the Ito Langevin equations of 
motion for system operators. In this picture, the state 
remains separable and the expectations will always have 



III we define these states and use them to the form of Eq. (16) and Eq. (18). 



derive master equations. 

II. FOCK STATE MASTER EQUATIONS 

In this section we derive master equations for a quan- 
tum system interacting with a field prepared in a Fock 



At this point we must mention an important techni- 
cal issue. The composition rules for the quantum noise 
increments, expressed in Eq. ([5]), are generally modified 
for non- vacuum fields [421 164j . However, it is shown in 
Appendix |B 2[ that the Ito table for Fock states is identi- 
cal to that for vacuum. This allows the techniques from 
input-output theory to be extended to Fock states. 
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A. Fock-state master equations for the system 

Recall the first step towards deriving the vacuum mas- 
ter equation, Eq. Q, was taking the expectation of 
Eq. (IT]) with respect to vacuum, i.e. ]Eo,o[rf-'^]- Anal- 
ogously! to derive the Fock-state master equations we 
must take the asymmetric expectations, i.e. Eq. ( |16[ ) or 
Eq. ( [T8| . The only exphcit field operators in Eq. ^ 
are the quantum noise increments dBt and dAt. Conse- 
quently the action of the quantum noise increments on 
Fock states is needed: 



dBt\n^) =c^^^^e(^)l"-l^) 



dAt\n^) = dBly^^{t)\n - 1^ 



(19a) 
(19b) 



In Appendi x |B 1| we show how to derive these relations. 
Equations ( 19 1 show how "reference" Fock states of dif- 



ferent photon number couple through the quantum noise 
increments. 

We are now equipped to derive the Fock-state master 



equations. From Eq. (18 1, we take the partial trace over 
Fock states for an arbitrary system operator X (E) /field, 
whose equation of motion is given by Eq. ([T]) . Doing so 
yields the Hciscnberg master equations: 



dt 



„,„(i[ff,X])-f tn™,„(£t[L]x) (20) 
+ ^/m^{t)wm-lA[L\X]S) 

+ V^le(t)l'^m-l.„-l(5^^5 - X). 

To extract the Sc hrod inger-picture master equations, 
we make use of Eq. E™,„[X(0] = Trsy.[f?}„,„(i)^]- 

Then, using the cyclic property of the trace, we can write 
down the master equations for the system state: 



dt 



Qm.nit) — ^ i[H, Qm.ii] + ^^[L]Qm,n 
+ V"^C(*)['S'^'m-l,n,i'^] 

+ \/mn\£,{t)\'^{SQjn-l,n-l 



(21) 



Pm-l,n-l J 



This set of coupled differential equations is the main 
result of this section. The initial conditions for these 
equations are: the diagonal equations Qn^n should be ini- 
tialized with the initial system state Psys, while the off- 
diagonal equations should be initialized to zero. In order 
to calculate expectation values of system operators for an 
iV-photon Fock state one needs only the top-level density 
operator qn,n- However, extracting qn.n requires prop- 
agating all equations between and N to which it is 
coupled. We note some special cases of Eq. (21 1 have 
been derived previously in Refs. [551 SSI [SD] however lit- 
tle intuition or physical interpretation was given to these 
equations. 



The master equations in Eq. (21 1 require further ex- 
planation. The diagonal terms, Qn,m £^re valid state ma- 
trices describing the evolution of the system interacting 



with an n-photon Fock state for n G {0, . . . , A''}. For 
example, when A'^ = we recover the vacuum master 
equation: dpo,o = ~i[H-, Qofi]dt + C[L\QQ Qdt, which is the 
only closed- form equation in Eq. ( pTj ). For A^ > 1, the 
diagonal equations couple "downward" towards the vac- 
uum master equation via the off-diagonal equations Qm,n 
where m ^ n. These off-diagonal operators are non- 
Hermitian of trace-class zero 26J; consequently they are 
not valid state matrices but do satisfy Qm,n = Qnm- 

The fact that the equations couple downward means 
that we need only consider a finite set of equations, which 
can be integrated numerically and in some cases, analyt- 
ically. For a field in an A^-photon Fock state there are 
(A + 1)^ equations. From the symmetry Qn.m — pJn.n, 
the number of independent coupled equations reduces to 
1{N +1){N + 2). 

Finally, we comment on the physical interpretation of 
these equations. Absorption of a photon by the system 
significantly changes a field prepared in a Fock state, so 
its dynamics are non-Markovian [551 HH]- This neces- 
sitates propagating a set of coupled master equations. 
(In contrast, for coherent states photons can be removed 
while leaving the field state unchanged and a single mas- 
ter equation suffices.) Before the wavepacket has inter- 
acted with the system ^(i) is zero and only the top level 
equation qn^n contributes to the evolution of the sys- 
tem. In other words, the system evolves solely under the 



terms on the first line of Eq. (21), which describe evo- 



lution from an external Hamiltonian and decay due to 
coupling to the vacuum. When the wavepacket begins 
to interact with the system, ^(t) becomes nonzero and 
the other coupled equations contribute to the evolution 
of the system. Then, the information fiow propagates 
upwards from go.o to gN,N because the equations couple 
downwards. 

So far we have discussed the dynamics of the system 
before and during the interaction. The last physically 
important observation is related to the correlation be- 
tween the system and the outgoing field during and after 
the interaction. Consider the case where £^{t) is bimodal. 
When the temporal spacing between the peaks is much 
greater than the characteristic decay time of the system 
and since ^{t) is zero at these intermediate times, the co- 
herence between the first peak of the wavepacket and the 
system is lost before the second peak begins to interact. 
Thus only the top-level equation must be propagated at 
these times, and the only nonzero terms describe external 
Hamiltonian drive and decay into the vacuum. When the 
temporal spacing between the two peaks is on the order 
of the system decay time or shorter, then the initial tem- 
poral coherence between the peaks can affect the system. 



B. Output field quantities 

In addition to system observables, we may also be in- 
terested in features of the output field [5S]. Consider 
a field observable Y{t) with initial condition Yito) = 
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^3ys (X" Y. We insert the Ito Langevin equation of motion 
for Y into the asymmetric expectations. Using Eq. (18), 
i.e. the partial trace zz7,„_„(y(t)), the result is operator- 
valued Heisenberg master equations. We focus here on 
expectation values, ¥)m,n[Y{t)], which are found by trac- 
ing over the system as well, as in Eq. ( 16 ). For two quan- 
tities of interest - photon flux and field quadratures - we 
produce a set of coupled differential equations initialized 



in the same way as Eq. (21 ) 



1. Photon flux 

The photon fiux is given by dA^ , which counts the num- 
ber of photons in the field in the infinitessimal time in- 
crement t to t + dt \0, Sec. 11.3.1]. The rules of Ito 
calculus are used in Appendix |A 2| to give the equation 



of motion for the output photon flux A° , 

dA™' = L^Ldt + L^SdBt + S^LdBl + SdKf (22) 



Taking expectations over Fock states using Eq. ( 16 ) 
yields an equation for the mean photon flux. 



{t)]=W.raAL^L\ 



'^ie(t)pE,„_i,„_i[5^5]. 

(23) 



The solution to this equation ]E[A°"*(t)] gives the inte- 
grated mean photon number up to time t. 



2. Field quadratures 

A Hermitian field quadrature Zt measurable via homo- 
dyne detection is described by 



Zt = e"l'Bt + e-'"l'Bl 



(24) 



Following the same prescription, the equation of motion 
for the quadrature after the interaction is 

= e''l'{Ldt + SdBt) + e-"''{LUt + SUb}). (25) 



Taking expectations over Fock states using Eq. ( |16[ ) gives 
the mean homodyne current, 

^E„,,„ [Zr' m = E,„,„ [e^^L + e-^*Lt] (26) 



C. General input field states in the same 
wavepacket 

So far we have considered the case where the input 
field is a "pure" Fock state. These results can be gen- 
eralized to field states described by an arbitrary combi- 
nation (superposition and/or mixture) of Fock states in 



the same wavepacket. As the Fock states span the full 
Hilbert space, they form a basis for arbitrary states in 
the wavepacket ^{t), 



Pfiold 



m,n— 



|n^)(m^|. 



(27) 



The coefficients are constrained by the requirements of 
valid quantum states: pfioid > 0, Tr[pficid] = 1 and 

Afield = /'field 



When the input field is described by Eq. (27) the sys- 
tem state is 



f?total 



(0, 



(28) 



where Qm,n{^) ^'^^ solutions to the master equations. 
Generating the full, physical density operator for an ar- 
bitrary field requires combining the appropriate solutions 
from the hierarchy of coupled equations in Eq. (21) with 
associated weights c,„_„. The Heisenberg master equation 
is found in the same manner, 



TOtotal(0 — E ' 



(29) 



Finally, the expectation value of a system operator X is 
given by 

Etotai[X(i)] = Trsys+ficid [QLUt)x\ (30) 
-Ec™,„E„^„[A(i)]. (31) 



This technique also applies to the output field quanti- 
ties in Sec. IIIBI Note that the definition of the Hilbert- 
Schmidt inner product, Eq. (151, gives rise to the conju- 

but not in Eq. (30). 



gate coefficients in Eqs. (28 
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III. 



GENERAL N-PHOTON MASTER 
EQUATIONS 



In many experimental settings multiple photons are 
not created in Fock states. Fock states are a subset 
of more general A^-photon states, which have a definite 
number of photons but an arbitrary SDF Indeed, 
a quantum tomography protocol for characterizing the 
SDF was recently proposed |65] and implemented [66] . 
This motivates the derivation of master equations for 
such fields. 

In a single spatial and polarization mode, a general 
A^-photon state is 



1^. 



N 



dull ■ ■ ■ duiN "0(1^11 ■ • ■ 1 '-^n) 



(32) 



Again we assume quasi-monochromatic wavepackets such 
that '(/'(■) is a slowly- varying envelope with respect to the 
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carrier frequency. Then, in the time domain a general 
A^-photon state can be written as 



N 



dti . . . dtN i^{ti, . . . , In) 
xb\h)...b\tN)\0) 



(33) 



These states are not amenable to our analysis directly. 
Thankfully, a formalism for dealing with such iV-photon 
states has been developed [TTl E7] . 

To describe iV-photon states we make use of the oc- 
cupation number representation developed by Rohde et 
al. [U], which we review in Appendix [C] Using Eq. ( [CSj ), 
Eq. ( [Ml can be written in a basis of orthogonal Fock 
states, 



N 



E 

ii<...<ij\ 



A,; 



(34) 



where \nk^^) is a normalized Fock state described by 
Eq. ( |13[ ) with rife photons in basis function ^^(i). Count- 
ing the number of subscripts on A in Eq. ( 34 ) gives the 



total number of photons N, and the value of any sub- 
script ik reveals the basis function that photon is in. 

In order to derive the master equation, we must first 
write down the action of the quantum noise increments 



on Eq. (34) 



k 

dk^N) = dB\ ^a(i)l^w-i) 



(35) 
(36) 



Each master equations couples to a set of equations enu- 
merated by the indices {m, n, fc}, which run over the or- 
thogonal basis functions. The total number of equations 
required to describe such a state depends on the overlap 
of the initial wavepacket with the particular choice of ba- 
sis. Equations for the output field can also be derived for 
iV-photon states, but we omit them for brevity. Equa- 
tions similar to Eq. ( 38 ) were derived in Ref. [55] for two 



photons but did not include dh.t or S. 

Finally, we can consider input fields in combinations 
(superpositions and/or mixtures) of different iV-photon 
states. In particular we allow the total state to be a com- 
bination of different states with the same photon number 
and a combination of states with different photon num- 
bers. To describe such a state first we need to consider a 
general combination of iV-photon states. That is. 



cl,k\Kn){Ln\ 

KXG{ip,4,..a} 



(39) 



where the summation is over different states with the 
same photon number iV. Then we can sum over photon 
numbers to obtain the most general input field: 



Pficld 



p=0 



(40) 



The coefficients Cl,k and Cp are constrained by the re- 
quirement that the input state be a valid quantum state. 
Using Eqs. (40 1 and (28), the equations for the system 



and output field can be found. 



where IV'at-i) is defined as 



k \ 

N-l) 



(37) 



IV. EXAMPLE: FOCK-STATE MASTER 
EQUATIONS FOR A TWO-LEVEL ATOM 
INTERACTING WITH A GAUSSIAN 
WAVEPACKET 



and is interpreted to mean that a single photon in one of 
the basis Fock states has been annihilated. 

To derive the master equation for a system inter- 
acting with the field \iPn), an asymmetric expectation 
value needs to be defined for such states: E^^_^^[0] — 
Tr[(yOsys ^ \4'n) {4'm\Y O] . As before this defines the gener- 
alized density operators Q^p^^ip^- Using these definitions 
the master equations for the generalized density opera- 
tors are 



dt 



(38) 



k 
k 

Y^/mknk'Ck{t)£.k'{t) 

k,k' 



Efficient photon absorption is important for informa- 
tion transfer from a flying to a stationary qubit. In this 
section we analyze this problem with a study of the ex- 
citation probability and output field quantities for Fock 
states interacting with a two-level atom. This problem 
has been studied before in much detail for a single photon 
in Refs. [36U38j . Our intention is to make a direct con- 
nection to established results and then to extend those 
results to higher photon numbers. Consequently, we do 
not focus on optimizing wavepacket shapes as other stud- 
ies have [36H39] . 

The single-mode approximation in Sec. [T| is rooted in 
the presumption that the wavepacket can be efficiently 
coupled to the two-level atom. This has been considered 
in the case of a mode-matched wavepacket covering the 
entirety of the An solid angle in free space [36^.37 . A more 
widely applicable context is that of strongly confined ID 
photonic waveguides [39]. In such systems the coupling 
rate into the guided modes T g can be much larger than 
into all other modes Fj^, where the total spontaneous 



emission rate is F 



In the following 



7 



analysis, we take the idealized limit that coupling to all 
other modes can be fully suppressed and we set T± — 0. 
To properly account for losses, a second mode can be 
introduced using the tools of Sec. |V] and finally traced 
over. 

In Sec. |IV A[ we examine the form of the master equa- 
tion for the simple case of a two-photon Fock state. Next 
in Sec. |IVB| we numerically examine a two-level atom 
interacting via a dipole Hamiltonian with a wavepacket 
prepared with at most two photons. First we reproduce 
the single-photon excitation results from prior studies, 
then we broaden these results to include two photons and 
output field quantities. Finally in Sec. IV C we present 



a numerical study for large-photon-number Fock states. 
This allows us to explore the relationship between exci- 



tation probability, bandwidth, interaction time, and pho- 
ton number. For photon numbers iV ^ 1, we identify a 
region of strong coupling. 



A. Two-photon Fock state master equations 

It is instructive to examine the form of the master 
equation for the simple case of interaction with a two- 
photon Fock state where both photons are created in the 
same temporal wave packet S,(t), |'0)fleid = |2{). From 
Eq. (21), the two-photon Fock state master equations 



Q2At) 


= C[L]g2,2 


- i[H, Q2,2] 


+ V2m[SQi.2.L^] + ^e{t)[L, Q2.iS^] + 2|e(t)P {Sqi^iS^ - gi,i) 


(41a) 




= C[L]g2,i 


- i[H, Q2,l\ 


+ V2m[SQiA.L^+C{t)[L, g2,oS^ + V2\m? {Sg^^S^ - q^^^) 


(41b) 






" Q2fi] 




(41c) 


QiAt) 




- i[H, Qi,i\ 




(41d) 




= i^[L]Qi,o 


- i[H, Qifi] 




(41e) 






- Qofi] 




(41f) 



with the initial conditions: 



Q2.M ^ Ql.M ^ QO,m ^ Psys (42) 

g2,l(o)-^^2,o(0)-^^l,o(0)-o. (43) 



Similar equations to Eqs. (41 1 were originally derived 
in Ref. [551 Equations 71 (a)^(f)] for a two-level atom 
but without the 5* operator and the term proportional to 
|^(i)p. For an arbitrary quantum system and single pho- 
ton equations which include S and the term proportional 
to |^(t)|2 were later derived in Ref. |49j. Then Ref. jSU] 
showed how to propagate these equations for any super- 
position or mixture of one photon and vacuum. 

Now suppose the input field is in a superposition of one 
and two photons, |V')fioid = ct\^i)+P\'^d with + — 
1. From Eq. (28) we combine the solutions to the master 



equations, Eq. (41), to get the physical state. 



^?total(^) 



\a\'Q2At) + WQi.At) (44) 
+ a*PQ2,i{t) + aP*Qi.2{t)- 



Notice that the last two terms of Eq. (44) originate in 
the coherences of the input field. It is interesting that 
the "off-diagonal," traceless, generalized density opera- 
tors (e.g. Qi^2) contribute to the calculation of physical 
quantities, albeit in Hermitian combinations. Had the 
field been a "pure" Fock state or a statistical mixture of 
one and two photons, these terms would not appear. 



Output field quantities are calculated in the same fash- 
ion as Eq. (|44|). For example, the mean photon flux is. 



lEtotai[Ar*(i)] = |arE2,2[A° 
-ha/3*E2,i[A 



out] 



out] 

t J 



a*/3Ei,2[Ar*], (45) 



where Eq. (30) was used to calculate IEtotai[-] 



A two-level atom interacting with one- 
two-photon Gaussian wavepackets 



and 



Now we specialize to a wavepacket prepared with up 
to two photons interacting on a dipole transition with a 
two- level atom initially in the ground state \g). In the 
absence of an external system Hamiltonian the master 
equation parameters are: H — 0, L — \/T\g){e\, S — I, 
and the coupling rate is chosen for simplicity to be F = 1. 
We focus on a square-normalized Gaussian wavepacket, 
as defined in Ref. [37], whose peak arrives at time ta, 



2^ 



1/4 



exp 



--rit-taY 



(46) 



with no detuning and frequency bandwidth f2. For Gaus- 
sian wavepackets the simple relationship between band- 
width and temporal width enables us to explore the 
tradeoff between interaction time and spectral support 
around resonance 1861. 
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time (tr) 



FIG. 2: Comparison of a Gaussian wavepacket of bandwidth Q,/T — 1.46 in three initial field states: a single-photon Fock state 
(solid), a two-photon Fock state (dashed), and an equal superposition (dash-dot). The wavepacket |5(f)P is shown in black 
filled grey, (a) Excitation probability of a two-level atom, (b) Photon flux. It is distinctly modified by interaction with the 
atom, (c) Integrated photon flux. For comparison the integrated single-photon flux is plotted when there is no atom. 



To study the excitation probability we numerically in- 
tegrate the master equations (41a)--(41f ). Then, for a 



given input field state we calculate the excitation proba- 
bility. 



Pe(0 -Tr[etotaiW|e)(e|], 



(47) 



where ptotai is given by Eq. ( 28 ) 



Figure ^a) presents the excitation probability for a 
two-level atom interacting with a Gaussian wavepacket 
Eq. ( 46 ) prepared in a "pure" Fock state of one and two 
photons as well as an equal superposition; a — 13 — l/\/2 
in Eq. ( 44 ) . In the simulations we use a bandwidth known 



to be optimal for single-photon Gaussian wavepackets: 
f2/r = 1.46 [36 . This gives a maximum excitation prob- 
ability of P™'*'' « 0.801 for TV = 1 as found in other works 
[551 - I55] . Putting a second photon in the wavepacket 
slightly increases this to p™^^ « 0.805; however, we see 
in Sec. IIV CI that this is not universal behavior for all 
bandwidths and photon numbers. 

In Fig.[2lb) we plot the mean photon flux of the output 
field, (iIE[/\("*]/(it, after interaction with the atom. For 
the single-photon wavepacket, we see a drastic change 
in the output photon flux when the photon is being ab- 
sorbed by the atom. For two photons, however, much 
of the wavepacket travels through the atom undisturbed, 
since a two-level atom can absorb at most one photon. 
The related integrated mean photon flux, IE[A°^'], is plot- 
ted in Fig. ^c). For these "pure" one- and two-photon 
Fock states there exist a definite number of excitations. 
Any excitation induced in the atom through absorption 
of a photon eventually decays back into the field. This 
is shown in Fig. where the integrated mean photon 
flux for long times approaches the number of initial ex- 
citations {1, 1.5,2}. During the absorption of the single- 
photon wavepacket, the integrated intensity flattens out 
since the photon has been transferred to an atomic exci- 
tation and arrives only later after decay. 



For a single-photon wavepacket, the Schrodinger equa- 
tion can be solved analytically for the excitation proba- 
bility [351138]: 



dt'm 



(48) 



The simulations in Fig. [2] agree with the analytic expres- 
sion in Eq. (48 1. However, it is not clear that the method 



used to derive Eq. ( |48[ ) can be extended to higher photon 
numbers. 



C. Excitation for large photon numbers 

In this section we expand the numerical study of exci- 
tation probability to Gaussian wavepackets of the form 
of Eq. ( 46 1 prepared Fock states with photon number 
iV > 1. 



1. Scaling 

For small bandwidths (f^/F ^ 1), see the left side of 
Fig. [sf^a) , one would expect a high probability of excita- 
tion from the substantial spectral support near the tran- 
sition frequency of the atom. However, the long temporal 
extent of the wavepacket means the photon density over 
the relevant interaction time scale r = l/F is too small to 
significantly excite the atom 37] ■ A complementary way 
of understanding this is that the dissipative terms in the 
master equations [terms on the first line of Eq. (21 1] pre- 
vail over the coherent coupling (terms on the other lines) . 
By extending the analysis in Ref. [26], we find a recur- 
sive scaling of the excitation probability for very wide 
wavepackets: PJ°^'' « Pn, where Pn = NPi{l - 2Pn-i) 
with Pi = 4max|^(i)p. 
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photon number (N) photon number (N) 

FIG. 3: (a) Maximum excitation probability P™'"' of a two- 
level atom interacting with Gaussian wavepackets of band- 
width n/r for photon numbers A'^ £ {1, . . . 10}. Small (large) 
bandwidths correspond to long (short) temporal wavepackets. 
(b) Scahng of P™'"' with photon number (red circles). The fit 
shown is Pf='''(Af) = 1-^0.269 iV"°-^^=' (blue line), (c ) Scaling 
of p™^'' with optimal bandwidth for each photon number A*' 
(red circles). The fit shown is nopt(Af)/r = IA5N°-^^'' (blue 
line). Details of the fits can be found in the main text. 



In the other asymptotic regime where bandwidths are 
large (r2/r ^ 1), see the right side of Fig. |3][a), the 
maximum excitation probability is small even for large 
photon numbers. This is due to the wavepacket being 
so short that its bandwidth is spread over frequencies 
far from the atomic resonance. We numerically find the 
asymptotic scahng P^'''' = 5NT/n for f^/F e [10^, 10^] 
with = 1 for photon numbers N G {1, . . . , 10}. 

At intermediate bandwidths, we note several interest- 
ing features. First, the maximum excitation probabilities 
are not universally ordered by photon number and adding 
photons to the field can decrease P"'*''. In fact, there ex- 
ists a bandwidth region in Fig. [3] where a single photon in 
the wavepacket is optimal for excitation, il/F « [.5, 1.4]. 

Second, for each photon number there exists an op- 
timal bandwidth for excitation. In Fig. |3] (b) we have 
plotted the absolute maximum of Pg (maximized over 
t and r2/F) as a function of the number of photons. 
We find excellent agreement (i?^ = 1) by fitting to 
the model IPf'^'iN) = 1 - aN^'' over the range N e 
{10, ...,40} with coefficients (95% confidence): a — 
0.2694(0.2678,0.271), b = 0.973(0.9709,0.975). There- 
fore the absolute maximum of Pe does monotonically in- 
crease with N, but with diminishing returns. 

In Fig. [3] (c) we investigate the optimal bandwidth for 
excitation for each photon number N. Fitting to the 
model r2'"^'^(iV)/F = aN'' gives a = 1.447(1.418,1.476) 



and b = 0.9869(0.981,0.9928) with 95% confidence and 
— 0.9998. Thus, to achieve this scaling for photon 
number N, the optimal bandwidth of the wavepacket is 
floptiN)/r « 1.45iV°-^s'^. Thus, the optimal width seems 
to be proportional to the single-photon optimal band- 
width, r2opt(A^)/r « 1.46iV. 



2. Dynamics 

Finally we illustrate the excitation probability dy- 
namics. Figure ^ shows Pg for bandwidths fl/T e 
{50,1,1/20}, chosen to illustrate three types of behav- 
ior. In each subplot (a)-(c), excitation curves are plotted 
for photon numbers N E {1, . . . , 10}. 

In Fig. [ijja) a short pulse quickly excites the atom, 
which then decays into vacuum with rate F after the 
wavepacket leaves the interaction region. Larger photon 
number corresponds directly to larger maximum excita- 
tion. In the intermediate bandwidth regime, fi/F « 1, 
excitations can be coherently exchanged between the 
atom and field, leading to oscillations in the excitation 
probabilities. This continues until the wavepacket leaves 
the interaction region as shown in Fig. |4]jb). Similar 
damped Rabi oscillations were observed for large-photon- 
number coherent state wavepackets in Ref. [371 Fig- 5]- 
For a single photon in the field, these oscillations are 
never seen due to the tradeoff between spectral band- 
width and photon density [27l |68] • At the chosen band- 
width Q/T — 1, a single photon achieves the highest 
maximum excitation with maximum excitation falling off 
roughly with photon number in agreement with Fig. [3] 
Finally, in Fig. |4]jc) we see that an atom interacting with 
a long wavepacket is excited and then decays well within 
the wavepacket envelope and the Pe(i) curves are nearly 
symmetric around the peak of the wavepacket for all pho- 
ton numbers = {1, . . . , 10}. 



3. Strong coupling 

The damped Rabi oscillations seen in Fig.|4]jb) suggest 
that there is a regime where coherent processes dominate 
over dissipation, known in cavity QED as the strong cou- 
pling regime. The authors of Ref. 68J defined a strong 
coupling parameter (for very short rectangular wavepack- 
ets): VNgoff ^ r where g^ff — ^(t) y^F^. Specifically the 
wavepacket was taken to be ^{t) = l/V^max for times 
t < ^max ^ l/F and zero otherwise. In this limit they 
showed that full Rabi oscillations for N photons occur 
at frequency uju = gcsV^- In Fig. [5] we compare their 
analytically-predicted excitation oscillations with our nu- 
merical calculations for = 50 photons. In (a), the 
wavepacket is long compared to 1/F and, while the os- 
cillation frequencies match, the amplitudes do not due 
to dissipation. For short wavepackets, as seen in (b), co- 
herent coupling prevails over dissipation, we see excellent 
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time [tr) 



FIG. 4: Excitation probability Pe of a two- level atom interacting with Gaussian wavepackets of bandwidth fi/F — {50,1,1/20} 
prepared with iV G {1,...,10} photons. Highlighted are A'' = 1 (solid), N = 2 (dashed), and iV = 10 (dash-dot). The 
wavepacket |C(i)P is plotted in black filled grey (normalized in (a) for clarity), (a) Behavior of short temporal wavepackets 
(large bandwidths) shows Pe is ordered by photon number, (b) For intermediate bandwidths, we see damped Rabi oscillations, 
discussed in Sec. |IV (Ts] Note that Pe is not necessarily ordered, (c) Behavior of long temporal wavepackets (small bandwidths) 
where Pe is again ordered. Note the different time scales in (a), (b), and (c). 



agreement with the predicted frequency (in our param- 
eters: ujR — 2£^{t)^/TgN) and good agreement with the 
predicted amplitude. 

For non-rectangular pulses the frequency of the Rabi 
oscillations is time-dependent as seen in Fig. gb). We 
must account for the time variation of the wavepacket 
^(t) in order to define a more general strong coupling pa- 
rameter. To achieve strong coupling, the coherent cou- 
pUng rate into the guided modes \/ NT g\S^{t) \ must dom- 
inate the total relaxation rate F. We can immediately 
define the condition for instantaneous strong coupling: 
yjNTg\£^{t)\/T ^ 1. However, in order to see interesting 
dynamics such as a complete Rabi oscillation, the cou- 
pling must remain strong over a characteristic timescale 
T. From this argument we define an average strong cou- 



pling parameter, 




5 . , , 0.04 0.08 0.12 
time (tr) 



FIG. 5: Comparison of the numerically-calculated (dark blue) 
and analytically-predicted (dashed orange) Rabi oscillations 
for rectangular wavepackets (normalized for clarity) with A'' = 
50 photons, (a) Wavepacket length tmax large compared to 
l/F. (b) Wavepacket length approaching the limit tmax ^ 
l/F. We see increasing agreement between prediction and 
our numerics. 



^ " I dt\^{t)\ > 1 Vt, 



Ft 



(49) 



ts-T/2 



If, for any wavepacket £,{t), there is a value oitg such that 
Eq. ( 49 1 is much greater than one, then average strong 



coup 



ing has been achieved over the time window r. 
A natural choice for r is the characteristic decay time 
of the atom, l/F. In Fig. |6][a) we present a contour 
plot of the average strong coupling parameter for Gaus- 
sian wavepackets prepared in a single-photon Fock state 
[N — 1). Ideal couphng to the guided mode is assumed, 
Fg — F = 1. We see that, for any bandwidth, maxi- 
mum coupling occurs when the time window is centered 
at the Gaussian peak (indicated by the vertical, dashed 
white line) and that the strongest coupling is achieved 
for O/F = 4. Note that although the average strong cou- 
pling parameter for a single photon never exceeds one, 
for larger photon numbers the factor can lead to sig- 
nificant coupling. In Fig. [6|^b) the excitation probability 
dynamics are shown for an optimal bandwidth fi/F = 4 
wavepacket. We see the appearance of damped Rabi os- 
cillations when the wavepacket has iV = 50 photons that 
are completely absent when only a single photon is in 
the field. For comparison, a wavepacket of bandwidth 
0/F = 2 is shown in Fig. |6][c). Even at this bandwidth, 
damped Rabi oscillations appear for TV = 50 photons, 
albeit with reduced contrast and frequency. 



V. TWO-MODE FOCK STATE MASTER 
EQUATIONS 

In this section we derive the master equations for 
a system interacting with an arbitrary combination of 
continuous-mode Fock states in two modes (spatial or 
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FIG. 6: (a) Contour plot of the average strong coupling pa- 
rameter for a Gaussian wavepacket prepared with a single- 
photon as a function of center of the time window (ts) and 
bandwidth Q/T (where r — l/F). (b) and (c): Excitation 
probability of a two-level atom interacting with a wavepacket 
of bandwidths Q/F ^ i for (b) and Q/F ^ 2 for (c). Only 
A'' = 1 and A'^ = 50 photons are shown. The normalized 
wavepackets are shown in black filled grey. 



6ijl (see |6ni 



the ith field mode, and the scattering operator Sij 
strained by: SikS^^^. = 6ijl and Sl-Sk. 
Appendix A], [70l Sec. IV] and [71] and the references 
therein for more details on multi-mode QSDEs). Note 
that the subscript t on the multi-mode quantum noise 
increments has been dropped for notational compactness 
in favor of the mode labels The multi-mode quan- 

tum noise increments are defined, 



t+dt 



ds 6i(s), and dA^ 



t+dt 



dsbUs)b,{s). 



(51) 



The composition rules for these quantum noises incre- 
ments under Fock state expectation are 



dBidB^ = S^jdt, dBidAjk 



SijdBk, 



dK^-jdBl 



SjkdBj. 



(52) 



B. Two- mode Fock states 

We consider the case where photons in mode one are 
prepared in a temporal wavepacket ^{t) and those in 
mode two are in the wavepacket ri{t). The two-mode 
Fock state with N photons in mode one and Q photons 
in mode two is, 



1 



N 



BliO Blirj) |0;0) 



where the operators bI{-) are defined in Eq. (12). 



polarization). This generalization allows one to consider 
wavepackets scattering off of atoms or addressing multi- 
ple dipole transitions, for instance. The analysis for two 
modes is conceptually identical to but algebraically more 
complicated than the single-mode case. 



A. Multi-mode Ito Langevin equations 

The evolution of a system operator driven by multiple 
quantum noises is given by the multi-mode Ito Langevin 
equation, 

dX = {i[H, X]+^d^ [L^]X^ dt + [L\ , X]S,jdB^ 

i 

+ 4 [A, U]dB] + (SlxSkj - S,jX)dA,j . (50) 

where the modes are labeled by the subscripts {i,j,k} 
and repeated indices are summed. H is an external sys- 
tem Hamiltonian, the operator Li couples the system to 



Two-mode Fock-state master equations for the 
system 



Here we specialize the multi-mode equations, Eq. (50) 



and Eq. (51), to two modes by restricting the indices to 



run over the mode labels {1,2}. In Appendix |D] we show 
how do this calculation for any number of modes. We 
introduce notation for representing asymmetric expecta- 
tions over two-mode Fock states. 



IEm,n;p,9[A(t)] — Tr 



sys+ficld 



(Psys «> I ) (mj I (g) I ) (p^ I ) ^A (t ) 

(53) 



which also defines the two-mode generalized density op- 



erators Qm.n;p,q in analogy with Eq. (17). The reference 



field state is written as a tensor product where the la- 
bels {to, n\ refer to mode one and {p, q\ to mode two. 
The two-mode Heisenberg master equations is found by 
taking field expectations over the equation of motion 
Eq. (50). Thus, the action of the quantum noises on 
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two mode Fock states is needed: 

dBi\n^\q,^) = dt^/n^{t)\n - 1^: q^j) , 
dB2\n^;q.rj) = dty/qr]{t)\n^;q - 1,,), 
dAii|n^;q,,> = dBly/n(,{t)\n - l^;g^>, 
dAi2\n^; q^j) = dBl^i]{t)\n^;q - 1^). 



(54a) 
(54b) 

(54c) 

(54d) 



The action of dA2i and dA22 are defined in a similar way. 



We then obtain the Schrodinger-picture master equa- 
tions with Eq. ( 53 1 and the cyclic property of the trace, 



(55) 



+ \/ mn\£_{t)\'^ {^SiiQm-l,n-l;p,qS\i — Qm-l,n-l;p,qj + (ySi2Qm,n;p-l,q-lS\2 — Qm,n-,p-l,q-lj 



+ ^/mq £,{t)lf [t) S^2Q7yi-l,n-p,q-lSli + ^/TVp C{t)vit) Siig„i,n-l;p~l,qS, 



where as before there is an implied sum over repeated 
indices. The initial conditions are 



^?m,n;p,q(0) — Psys 
^?m,n;p,q(0) 



\i m — n and p ~ q (56) 
li m ^ n ov p ^ q. (57) 



To solve a two-mode master equation with N pho- 
tons in mode one and Q photons in mode two, pfieid — 
\N^){N(\ (g) we need to propagate {N + 1)^ x 

[Q + 1)^ coupled equations. As in the single-mode 
case the symmetries in the generalized density operators, 
reduce the number of independent 



0^ 

-7n,n;p,q ' 



equations to i(A^ 4- l)(iV + 2){Q + 1){Q + 2). 



D. General input field states in the same 
wavepacket 



system state is given by 

gtota\{t) = ^ C'^,n;p,q9m,n-p,q{t), (60) 
m,n,p,g 

where Qm,n\ p,q{t ) are the solutions to the master equa- 
tions in Eq. ( 55 ) . The composition for expectation values 
is given by 

]Et0tal[-'^(t)] = ^ C„i,„;p,gE [X{t)]. (61) 



As before, the conjugate coefficients in Eq. (601 come 
from the Hilber-Schmidt inner product, Eq. (15). This 



technique also applies to the output field quantities in 
Sec. Ell 



So far we have only considered the case where the input 
fields in mode one and two are in "pure" Fock states, 
although we allowed for different wavepackets. These 
results can be generalized to field states described by an 
arbitrary combination (superposition and/or mixture) of 
Fock states. Consider the state 

oo 

Pfiold = ^ Cjn,n-p,q\n^){me\®\qr,)(j)n\ (58) 
m,n,p,ij— 

oo 

= ^ Cm,n-p,q\n(^;qn){m^;Pn\. (59) 

m,n,p,g— 



E. Two-mode output field quantities 

The output field equations for two modes are signifi- 
cantly more complicated than the single-mode case be- 
cause one can consider linear combinations of the modes. 
Thus, there is a continuum of possible of output photon 
fluxes and fleld quadratures. Here we focus on photon 
flux and field quadrature observables that are diagonal 
in the modes. More complicated output observables that 
combine both modes can be obtained using beam splitter 
relations - effectively, a change of basis - as described in 
Ref. [7D]. 



As before, the coefficients, Cm,n-p,q, are constrained by the 
requirements of valid quantum states. For example the 
entangled NOON state for one photon is given by pscid = 
i(|le; 0)(l4;0| + |l4;0)(0; l^-f |0; l^)(le; 0| + |0^l4)(0; 1^1). 



When the input field is described by Eq. (58), the total 



1. Photon flux 



The number of photons scattered from mode j into 
mode i in the interval t to t + dt is given by (iA°"' . Its 
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equation of motion is 

dA°f = LlL.dt + ^S,,dB, + SlL.dBl + SlSi.dA,,. 

(62) 

Any possible two-mode photon counting distribution is 
given by taking expectations of Eq. (62). For example, 



With a straightforward generalization of the formalism 



developed in Appendix (|C|) and Sec. Ill one can derive 



master equations for states of the form of Eq. ( 66 ) 



Even more general is an A'^-photon state distributed 
over two modes bi and 62, 



tracing over the system and field for dAn gives the mean 
photon flux in mode one, 



1^ 



N 



duji . . . duN ^nI^^i, ■ ■ ■ ,^n) 



N 



d_ 
~dt 



ra,n;p,q 



(63) 



J{{a,h\{uj,)^P,hl{LOi))\Q). (67) 



where ai and j3i are weights for modes one and two, 
respectively. For example, if we set all the = in 



+ ^MC(^)IEm-l,n;p.g[i^<S'll] + ^/p7?(^)Em,r^;p- l,g [ij 'S'12] 

[S'jj-^i], Eq. (|67|) then there would be A^ photons in mode two. 

+ ^ \/rrm\£,{t)\^¥.rn-i,n-i-p,q[SliSji]. 



The equation for mode two follows similarly. 



2. Field quadratures 



The output quantum noise in mode i is given by 
= S,jdBj+L,dt. 



(64) 



Just as in the single-mode case, field quadratures are Her- 
mitian combinations of Bi and Bi^ . For instance, the 
field quadrature in mode one, = e^'^B^ +e~"^B/. The 
equation of motion for the mean ouput field quadrature 
Z™*, or homodyne current, after the interaction is. 



di 



E™,„;p,,[^r'(t)] =E,„,„;p,Je''^Li + e-'*Ll] (65) 



-f e"'"^^77*(<) E,n,n;p,q-l[Sli]. 



The equations for Zj"' follow similarly. 



F. General two-mode A'^-photon states 



The formalism developed in Sec. |III| suffices to describe 
arbitrary states in each mode separately and thus is di- 
rectly applicable to the two-mode master equations. 

A slightly more general case is when there are m pho- 
tons in mode one and N—m photons in mode two with an 
arbitrary spectral distribution function (such two-mode 
states can be entangled in the spectral degree of free- 
dom). These states can be written 



IV'at) = / ■■■J duji . . .dujN ^NiuJi, ■ ■ ■ ,ujn) (66) 
xb\{u}i)... b\{uJm)bl{uJm+i) . . . 4(wAr)|0) 



For a small number of photons it is tedious, but possi- 
ble, to write down the occupation number representation 
of the state in Eq. (67). Finding an efficient representa- 



tion for such state with arbitrary A^ is an open problem 
and would allow a derivation of general two-mode master 
equations. 



VI. TWO-MODE EXAMPLE: FOCK-STATE 
SCATTERING FROM A TWO-LEVEL ATOM 

In this section we illustrate the use of our two-mode 
formalism by examining the photon flux of the trans- 
mitted and reflected fields when one- and two-photon 
Fock states are incident on a two-level atom [27H3T| l33l - 
I35j . The two modes are the forward- and backward- 
propagating fields, as in a tightly-confined waveguide 
QED setting [3TJ [31] . As before we specialize to a Gaus- 
sian wavepacket ^(t) described by Eq. (46). The master 
equation parameters we use are again those for dipole 
coupling without external Hamiltonian drive: H ~ 0, 
Li = A/Tils) (e|, Su = /, Sij = for z 7^ j, and the 
coupling rate is chosen to be = 1/2. The forward- 
propagating field is prepared in a Fock state with N E 
{1, . . . , 5} photons while the backward mode is initially 
in vacuum; that is, |V-'ficid) = |A^c;0)- 

In Fig. [7](a) we plot the excitation probability Pg for a 
two-level atom interacting with a wavepacket with band- 
width n/r = 1. The photon flux of the transmitted and 
reflected fields is plotted in Figures [tJc) and (d), normal- 
ized to the number of input photons A^. 

We first examine the single-photon input state (solid 
green curves). While absorbing the photon, the atom has 
a substantial Pg. The two peaks in the transmitted flux 
correspond to the attenuated input wavepacket and the 
contribution from remission into the forward mode 28]. 
Notice the dip between the peaks occurs when there is 
a large atomic excitation. Consequently this dip in the 
transmitted photon flux is due to atomic absorption and 
destructive interference with the incoming wavepacket 
[Si US [331 135]. Conversely, energy from the fleld that 
is not absorbed is scattered into the backward mode 
through the reemission process [55]. For A^ > 1, we see 
that the excitation probability is comparable to that for a 
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single photon, but the relative transmitted and reflected 
photon fluxes are quite difi^erent. In particular the ratio 
of transmitted to reflected flux increases with N. 

In order to understand this phenomena it is necessary 
to consider the normalized transmitted and reflected pho- 
ton numbers in the long-time limit (]E[Aii]and IE[A22]) 
at different bandwidths |28l [34] . In Fig. ^ we explore 
this issue numerically. Recall that the reflection process 
is facilitated by absorption and then reemission into the 
backward mode. Thus one would expect reflection to 
dominate for small bandwidth wavepackets, which is in- 
deed what is seen in the left hand side of Fig. [8j In 
the large bandwidth limit very little of the wavepacket 
is near resonance with the atomic transition so no ab- 
sorption occurs and the wavepacket is transmitted. The 
bump in the > 1 transmission and reflection curves is 
a consequence of an effective photon-photon interaction, 
explained in Refs. [SU |34] . By calculating the scatter- 
ing eigenstates, Zheng et al. found "multi-photon bound 
states" |34| which can increase transmission in that band- 
width region. 

It is also possible to examine scattering between the 
forward and backward modes, as was studied in Ref. [34] . 
by propagating the equations for A12 and A21; however, 
we omit this analysis for brevity. 
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FIG. 7: Scattering of a Gaussian wavepacket of bandwidth 
ft/T = 1 from a two-level atom. The wavepacket (black 
filled grey) is prepared with G {1,...,5} photons, (a) 
Excitation probability. Photon flux of the transmitted (b) 
and reflected (c) fields, normalized to input photon number. 
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FIG. 8: Normalized transmission and reflection for Gaus- 
sian wavepackets, prepared with N € {1, • • • , 5} photons, with 
bandwidths Q./V scattering from a two-level atom. The left 
(right) side represents long (short) temporal wavepackets. For 
larger photon number, note the increased transmission at in- 
termediate bandwidths. 



VII. DISCUSSION 

In this paper we have derived master equations for an 
arbitrary quantum system interacting with a continuous- 
mode Fock state in one or two modes (spatial and/or 
polarization). We generalized these results to include 
superpositions and mixtures of A^-photon states with ar- 
bitrary spectral distribution functions, and thus we can 
describe interaction with very general states of light. 

The power of our formalism lies in its direct applicabil- 
ity to more general systems of interest in quantum optics 
such as multi-level atoms, symmetrically-coupled atomic 
ensembles, and continuous variable systems such as nano 
mechanical resonators. For example, it is possible to re- 
produce the cavity-mediated, single-photon pulse shap- 
ing results of Ref. First we identify that H — 0, 
L v'ra, and S — I are the relevant substititutions. 



Then, our expression for the output photon flux, Eq. ( 23 1 , 
is equivalent to Equation (22) in Ref. [63] for one photon 
(i.e. in our equations set A^max = 1) after some algebraic 
gymnastics. 

As pedagogical examples, we studied features of Fock 
states interacting with a two-level atom in one and two 
modes. In the single-mode model [Sec. IV we saw the 
maximum excitation probability P™^'^ was low for both 
smaU {n/r < 10-1) and large [n/T > 10^) band- 
widths. The low p™^'^ for small bandwidths, centered at 
the atomic resonance, might seem counter intuitive. In 
the time domain the corresponding wavepacket is broad, 
nevertheless the near-resonant photons all get absorbed, 
but are immediately reemitted by the vacuum coupling, 
which leads to a small average Pg. This intuition is con- 
firmed in the two- mode simulations, presented in Scc. |VI[ 
where wavepackets with small bandwidths are nearly per- 
fectly reflected. The reflection is mediated by photon 
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absorption and the consequent reemission, which is di- 
rectionally unbiased. However, destructive interference 
between the incoming wavepacket and the transmitted 
mode results in reflection only; i.e. the atom can act as 
a perfect reflector. 

A detailed investigation of this phenomenon requires 
access to the individual quantum trajectories [75] rather 
than the ensemble averaged evolution given by the mas- 
ter equations. The differential equations for the quan- 
tum trajectories, known as stochastic master equations 
or "quantum filters" |73j . for a single photon in one 
mode have been derived previously for photon-counting 
[26t l49| [50] and homodyne [49l [50] measurements of the 
output fields. We are presently extending these to Fock 
states in one and two modes. Access to the conditional 
states would allow for measurement-based feedback con- 
trol [M] , 

A number of interesting applications of our formal- 
ism remain to be explored, including the investigation of 
pulse shaping for few-photon states, high efficiency quan- 
tum memories, and mediated photon-photon interac- 
tions. Our formalism is particularly applicable to quan- 
tum networks [H [5]. Recently, the theory of cascaded 
quantum systems [JT] HH] has been formalized to the 
point where simple rules for composing modular quan- 
tum optical systems into a network have been developed 
[Sll [TO] [m HH [75] . One needs only the (S", L, iJ )-tuple of 
each module specified in order to perform network anal- 
ysis and simplification. As our description of the system, 
input, and output fields is also in terms of a (S, L, H)- 
tuple, it is likely that our formalism can be ported to this 
setting. 
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Appendix A: Quantum noise and quantum 
stochastic calculus 



A rich mathematical machinery forms the foundation 
for the manipulation of QSDEs and their derivation from 
physical systems. Here we only touch the surface com- 
mensurate with our purposes; an interested reader is di- 



rected to Refs. [m [H [M [351 [57] [Sa [Ml [ZBHZS] for a 

more rigorous and detailed analysis. 

We present an introduction to the formalism of quan- 
tum stochastic calculus through the canonical example 
of a two-level atom interacting with a quantized, one- 
dimensional field. The atomic raising and lowering oper- 
ators are tT_|_ = \e){g\ and cr_ = \g){e\ with transition fre- 
quency Wo . The field is described by creation and annihi- 
lation operators, a\uj) and a{uj), obeying the commuta- 
tion relation [0(0;), a'^ (w')] = i5(a; — w'). The interaction- 
picture coupling between the atom and the field, within 
the rotating wave approximation, is 



dLu K{uj)aiu)e~''^'^"^'>^' + II.C., (Al) 



where the dipole coupling, k{uj) = \{e\d\g)\y^uj/4:'KeoficA, 
has units of \/ frequency and A is the effective transverse 
cross-sectional area of the mode, see Domokos et al. 1271. 



1. The quantum white noise limit 

To take the quantum white noise limit wc first assume 
weak coupling, i.e. that |«;(tLi)p ^ Wq. When k{uj) is 
slowly varying around loq, we make the Markov approx- 
imation that the atom has a flat spectral response; i.e. 
k{lo) — k{ujq). This implies that the correlation time 
of the field is short compared to the slowly-varying in- 
teraction time, Ts ~ l/|K(ajo)p. From the perspective 
of atomic operators, the field is (5-correlated in time and 
retains no memory of its past interactions. In this limit 
we can introduce the following field operators. 



b{t) = 



1 



(A2) 



which obey the commutation relation [b{t), b^{t')] = 
6{t — t'). For classical stochastic processes, 5-correlation 
implies white noise, so the operators h{t) and h^{t) are 
dubbed quantum white noise operators. Recast in terms 
of these operators the interaction Hamiltonian is 



HintW = i^{a^b\t)-a+h{t)) 
where we define k{ijJq) = \/ ^j^i: and set h—\. 



(A3) 



Under the white noise-driven Hamiltonian in Eq. ( A3 1 , 



the system and the field undergo joint unitary evolution 
via the propagator U{t^ that satisfies the Schrodinger 
equation. 



dt 



U{t)^^{o^ bHt)-a+b{t))U{t). 



(A4) 



This expression defies rigorous mathematical definition 
due to the singular commutation relation of the opera- 
tors b{t) and b^{t). To remedy this we first consider the 
quantum stochastic processes. 



ds b{s) and Bj 



ds b''{£ 



(A5) 
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The singular nature of the quantum white noise oper- 
ators can be removed by expressing Eq. ( A4 1 in terms 



of the continuous differential increments dBt and of 



Eq. (A5): 



/t-\-dt nz^az 
dsh{s)^dBt and dsb\s)^ dB\. (A6) 

These are the quantum, non-commuting analogues of the 
classical Wiener process and are referred to generically 
as quantum noise increments. Now equation ( A4 ) can be 
recast in differential form: 



t+dt 



dUt = V7 (cj- dBj - CJ+ dBt) o [/, 



(A7) 



Although technically an integral equation, this is referred 
to as a quantum stochastic differential equation (QSDE). 

In contrast to ordinary differential equations, white 
noise QSDEs have equivalent but non-identical represen- 
tations. Equation (A7| is an example of a Stratonovich 



QSDE, identified by the notation dBtoUt, which indicates 
the ordering of dBt and Ut is important; i.e. that they do 
not commute. Stratonovich QSDEs arise as the natural 
form for the quantum white noise limit of physical pro- 
cesses [7B] and follow the rules of standard calculus. More 
amenable for our purposes is the Ito form of a white noise 
QSDE. The quantum Ito integral is defined such that the 
integrand and the operator differential, dBt, act on inde- 
pendent time intervals and therefore commute, which is 
useful for taking expectations. Thus, we will work exclu- 
sively with QSDEs in Ito form, denoted simply by dBtUt- 
However, the Ito form brings the burden of its own calcu- 
lus, which requires that differentials be taken to second 
order. 

Performing the conversion from Stratonovich to Ito 
form [421 176] on Eq. (A7| and renormalizing a trivial 
energy shift, we obtain the QSDE for the unitary time- 
evolution operator 



dUt 



iy^a^ dBj - y/ja+ dBt - dt) Ut- 



(A8) 

The first two terms represent the atomic dipole coupling 
to the quantum noise increments, and the third deter- 
ministic term is an artifact of the transformation from 
Stratonovich to Ito form, known as the Ito correction. 



2. General stochastic time evolution operator 

The quantum white noise limit can be extended to in- 
clude coupling of a system operator S to the number of 
photons in the field at time t. This interaction Hamilto- 
nian is 



H^nnM=~Sb\t)b{t). 



(A9) 



From this Hamiltonian we identify a third fundamental 
quantum noise which can drive the system in the white 



noise limit, 

At- 

which has increments 

(•t+dt 



ds b\s)h{i 



dsb^s)b{s) dAt 



(AlO) 



(All) 



Including the possibility of an external system Hamilto- 
nian H, the most general QSDE for the time evolution 
operator in one mode has the form |70| 

dUt = {^L^L + iH)dt- L'fSdBt (A12) 

+ LdBl + {S~I)dAt}Ut. 

This equation describes the coupling of system operators 
L, L\ and S to the quantum noises dB\, dBt, dAt, and 
/ is the identity operator. The system operator S can 
be found from the bare Hamiltonian coupling of S in 
Eq. [A9\ with rules described in Ref. [SD|. 



3. Ito Langevin equations 



The time evolution operator in Eq. (A12| allows us to 
calculate the equation of motion for an operator O. Since 
we work with Ito QSDEs, this requires taking differentials 
to second order, 

d{u}OUt) = dUlOUt + ulOdUt + dU^OdUt- (A13) 

Note that in the literature one may encounter the "quan- 
tum flow" notation where an operator O at time t is given 
in the Heisenberg picture by jt{0) = UjOUt- When 



manipulating QSDEs such as Eq. (A13) one encounters 
products of the quantum noise increments. Under vac- 
uum expectation the rules for these products are given 
by the vacuum Ito table 



X 


dBt 


dAt 


dBl 


dt 


dBt 





dBt 


dt 





dAt 





dAt 


dBt 





dBt 














dt 















(A14) 



where we take the row and multiply by the column (row 
X column) to obtain the resulting product under vacuum. 
With Eq. (|A13| and Eq. ^AT^ we can write down the 



Ito QSDE for an operator X ® /fit 
dX ={i[H, X]+ C''[L]X)dt + 



Id, 



[L\X]SdBt 
S''[X,L]dBl + (S^XS - X)dAt, 



(A15) 



referred to as an Ito Langevin equation. Further, we can 
write down the Ito Langevin equation for output field 
quantities, such as the quantum noise 



Ldt + SdBt 



(A16) 
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and photon number A™ 



2. Fock and A^-photon Ito tables 



dA° 



L^Ldt + S^LSdBj + L'^SdBt + S'^ SdKf 



(A17) 



4. Multi-mode time evolution operator 

The evolution of a system driven by muhiple quantum 
noises is given by the QSDE for the multi-mode time 
evolution operator, 

dUt ={(Sy - %/)dA,j - L\S,,dB, + UdB\ 

-{\L\L,+iH)dt^Uu (A18) 

where is the coupling between the ith mode and the 
system, H is an external Hamiltonian, and the scatter- 
ing operator Sij is constrained by: SikSj/, = SijI and 

Sl^Skj = S^jl (see [31 Appendix A] and [TDl Sec. IV] 
and the references therein for more details on multi-mode 
QSDEs). Note that the subscript t on the quantum noises 
has been dropped for notational compactness in favor of 
the mode labels {«, j}- The multi-mode quantum noise 
increments are defined: 

/t+dt I't+dt 
ds bi{s) 1-^ dBi, and / ds 6|(s)6j(s) M- dA^. 

(A19) 



Appendix B: Quantum stochastic calculus for Fock 
states 

1. Action of the quantum noise increments on Fock 

states 

Recall the single photon state is defined by |1{) = 
/ ds ^(s)&t(s)|0) = B\()\0). Acting the quantum noise 
increment dBt on this state gives 

/t + dt n 
drb{r) J ds^sPisM 

t+dt 



drds {bHs)b{r) + 6{s - r))C(s)|0) 



pt+dt 



dtmio). 



(Bl) 



Some of this algebraic manipulation can be simplified by 
using the Gardiner-CoUett heuristic dBt = dtb{t) [i^ . 
Using this and the commutation relation [6(t), _B'''(^)] = 
^^(C)^(^) +^(^) this procedure is extended incrementally 
to higher photon numbers. Through induction we obtain, 

dBt\n^) ^ dty/^£,{t)\n - 1^>. (B2) 

By the same procedure we find the action of dAj, 

dAt\n^) = dB\^fnS,{t)\n - 1^). (B3) 



The vacuum Ito table, Eq. ( A14[ ), can require modi- 
fication for non-vacuum fields, such as thermal, coher- 
ent, and squeezed fields Here we show, sur- 
prisingly, that the Ito tables for continuous-mode Fock 
states and A'^-photon states are identical to the vacuum 
Ito table. This property was derived by the authors of 
Refs. [ini ISO] for a single photon, although never explic- 
itly written down in those papers [5T| . 

Consider the expectation of dBt dB\ for a single- 
photon Fock state. Normally ordering and simplifying 
gives 



m\dBtdBi\^) 



(1^ I {dBjdBt + dt)m) = dt. (B4) 



Alone, Eq. ( B4 ) is not enough to specify the Ito rule rule 



for for dBt dBj because the action of the noise incre- 
ments on Fock states couple different photon numbers, 
as in Eq. ( B2 ) . Consequently, we must consider cross 



expectations. Only after showing that {l^\dBtdBl\l^), 
{0\dBtdBlm), m\dBtdBl\0), and {0\dBtdBl\0) are pro- 
portional to dt can we say that dBt dBj — dt for the 
single photon Ito table. 

Now consider Fock states. One must show that 
{m^\dBtdBl\n^) = dt for all m and n. Thankfully it 
is straightforward to show that after normally ordering 
the operators - dBt,dB\,B{(^), and B'^ - the only sur- 
viving term is proportional to dt (terms proportional to 
dt^ are set to zero). Repeating this prescription for ev- 
ery product of the quantum noise increments in Eq. ( B4 1 , 



one can show the equivalence of the Fock and vacuum Ito 
tables. 

The Ito table for an A^-photon state with an arbi- 
trary spectral distribution function (within the quasi- 
monochromatic approximation) is also identical to the 
vacuum table. This follows from the occupation number 
representation, presented in Appendix [C] which relies on 
a decomposition in a basis of orthogonal Fock states, each 
of which respects the its own Fock Ito table. 



Appendix C: Occupation number representation for 
general A-photon states 

Here we review the occupation number representation 
of a general A^-photon state presented in Ref. [TT]. In 
one dimension and in a single mode, a general quasi- 
monochromatic A^-photon state can be written as 

IV'Jv) = duji . . . dujN ^j{oJi, ■ ■ ■ ,ujn) 



(CI) 

x6t(c^^)...5t(a,^)|0). 
In the time domain, this becomes 

IV'w) = J dti . . .dtN i^iti, . . . ,tN)b\ti) . . .b^t^m , 

(C2) 
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where the temporal envelope ^/'(ti, . . . , ^at) is the Fourier 
transform of ■ ■ ■ ,<jjjv) |12j ■ The temporal envelope 

is in general neither factorable nor symmetric in tk- It 
can be expanded in a set of complex- valued, orthonormal 
basis functions that satisfy J dt £,i — 



■0(^1, . . . ,iAr) 



E 

ii,...,ijv 



K,,...,^J^Ah).■.^^AtN). (C3) 



Each subscript runs over the labels for the basis func- 
tions, i.e. ik G {1, 2, . . . }. The expansion coefficients are 
given by the projection of the temporal envelope onto the 
basis functions, 



K,fs...x = J dh...dtNCih)...Cc{tNMh,...,tN). 

(C4) 

Defining a creation operator for a single photon in basis 
mode as B^(Ca) = / dt^a{tp{t), and using Eq. 

( C2|C4 |, we write the iV-photon state as 



\^n)= E K,,....^.BHC^J■■■BHC^,m. (C5) 
il,...,ijv 

Acting these operators on vacuum yields an expression 
for the A^-photon state in terms of basis Fock states, 
Eq. (13), in the basis functions. 



1^^) = 



A', 



.\/niln2l ■ 



'\n,^,)\n2^J... (C6) 



Counting the number of subscripts of A' gives the to- 
tal photon number N, which can be distributed among 
the basis Fock states in Eq. (C6). The number of pho- 
tons Ua in a particular basis function ^ait) is found by 
counting the number of indices of A' that are equal to a. 
For example, since they have 3 indices, the coefficients 
{AJ^ 12 describe a 3-photon state. The coefficient 
X[ I 4 refers to the state |2jJ|l^^), in which the first and 
second photons are in ^i(t) and the third in £,4{t). Due to 
the indistinguishability of photons, X'l 4 i and A4 ^ are 
also coefficients for the state 12^^)11^^), although they 
need not have the same value. In general, A^ is not 
invariant under permutation of its indices. The degree to 
which index-permutations are equal specifies the level of 
symmetry in the temporal envelope ip{ti, ...,tAr) [67ll82j. 
Following [TT], we define a new set of coefficients 



cr(ii,...,ijv) 



(C7) 



that sum over all permutations cr (in the symmetric group 
Sn) of the indices of coefficients of the type in Eq. ( C4 ) so 



that no two coefficients in Eq. ( C7 1 refer to the same basis 



Fock state. The iV-photon state of Eq. (C2|, written in 
terms of these coefficients, is 

IV-w) = E ^M,...,^JvhlCl)l"2^2)•■■ (C8) 

ii <-■ ■ <Z7v 

Now it is clear that these algebraic acrobatics have culmi- 
nated in a set of expansion coefficients that are precisely 
probability amplitudes. 



E |Aji,...,'ijv I" — 1: 



(C9) 



ii <---<Ziv 



and Eq. ( C8 ) is the occupation number representation of 



the general A^-photon state in Eq. ( C2 1 



Appendix D: Multi-mode expectations 

In this section we extend our formalism to a countable 
number of modes. First we define a multi-mode Fock 
state in T modes: 



I • a I • ■ ■ I ^ ' cj / 



1 



Vmi...NTl 



B\ur io), 

(Dl) 



where there are A^^ photons in the first mode with the 
envelope a{t) and ds |a(s)P — 1. 

To derive multi-mode mode master equations we must 
introduce notation, different from the main text, for rep- 
resenting asymmetric expectations. We define the multi- 
mode asymmetric expectation to be 



(D2) 



:Tr, 



sys4- field 



[psys®\ni){-m^o\®- ■ ■ ® 



ETrs 



where the superscripts and m} on E[.] (and g) refer 
to "reference states" in mode one. Note that Eq. (D2) 



also defines the generalized multi-mode density operators 

1 T 

n ;...;n /,\ 
; . .. ^ ' ' 

The final ingredient needed to derive the multi-mode 
mode master equation is the action of the quantum noise 
increments on Fock states: 



dA 



.nl) = dtVnieit)\ni;. 
.nl)=dBlV^9{t)\ni; 



1 



3 . 
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